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the coherent potential approximation. We discuss the T c as a function of No within 
the BCS, the McMillan, and the Belitz theory. The simplified Belitz theory gives the 
best agreement with experimental data. Since the density of states follows a simple 
power-law for accessible doping concetrations x, the present theory offers an analytical 
formula for T c {x). 
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1. Introduction 

A diamod doped with boron can become superconducting at temperatures of few kelvins 
PQ. Such critical temperatures are surprisingly high for impurity band conductivity. 
Moreover, one can expect that T c will be soon further increased. It would be useful 
to predict an optimal sample composition from a microscopic theory. To this end 
we evaluate the density of states at the Fermi level within the dynamical cluster 
approximation and show that it can be approximated by a simple power-law. Using 
this density of states the T c can be predicted from a McMillan type formula. 

Diamond is an insulator with a band gap of ~ 5.5 eV. Boron atoms create shallow 
acceptor levels close to the valence band with activation energy of holes ~ 0.37 eV. 
At a doping concentration above ~ 4.5 x 10 21 cm -3 the system presents a metal- 
insulator transition [2J and might become superconducting at 0.5 K. With increasing 
concentration T c increases. A favoured theory of superconductivity driving mechanism 
is, that electrons mostly interact with localized vibrational modes on the boron atoms 
in an optic spectra. However other possible theories are discussed as well [3]. 

Experimental data on boron doped diamond are controversial. Majority of samples 
were prepared in thin layers and their properties strongly depend on crystallographic 
orientation of the surface. There are two widely studied types with 100 and 111 
orientation, see references [2] and [I], respectively. Both groups have used the Microwave 
Plasma assisted Chemical Vapor Deposition method though under different growth 
conditions. Their samples differ in many characteristics including the T c . The 111 
samples have higher T c in general. In this paper we focus on the 100 samples which 
are also supported by bulk samples prepared with the High Temperature High Pressure 
method [1] revealing a comparable T c . 

The theory of boron-doped diamond can benefit from an extensive experience 
with superconductivity in disordered materials. [3 El [7J [8] Unlike in metallic alloys, 
where the disorder mainly modifies the coherence length and can be handled within 
simple approximations, in the case of doping induced superconductivity a sophisticated 
treatment of the disorder is necessary to describe a formation of the impurity band. 
Shirakawa et al [9] have shown that a concentration dependence of T c measured in 
boron doped diamond cannot be reproduced by the BCS theory in which the disorder 
is described on the level of the Coherent Potential Approximation (CPA). Since this 
theory yields too low values of T c , they proposed that disorder effects beyond the CPA 
are responsible for this disagreement. We will show that corrections beyond the CPA 
do increase T c , however in the extent which is not sufncent to cover discrepancies in 
question. 

In the present paper we compute the density of states at the Fermi energy level 
for boron doped diamond within the Dynamical Cluster Approximation (DCA). While 
the CPA assumes a single atom embeded in the effective medium, the DCA generalises 
this idea by embedding a cluster. The cluster allows us to describe splitting of closely 
located impurity bound states which contributes to the shape of the impurity band. 
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With respect to the superconductivity it is essential that the density of states at the 
Fermi level computed within the DCA is higher than the CPA value. 

Apparently, the density of states at Fermi level is much lower in the impurity 
band than in ordinary metals. Since values of T c are comparable, the boron doped 
diamond belongs to the family of the strong coupling materials. One can expect large 
dicrepancies from the BCS theory and we will show, that this is indeed the case. In this 
paper we use the Belitz theory [8] which generalises the McMillan formula to disordered 
sup er conduct ors . 

The structure of the paper is as follows. In Sec. II we simplify the Belitz theory 
so that it depends on the doping exclusively via the disorder dependent density of 
states. In Sec. Ill we evaluate the density of states within the dynamical cluster 
approximation. The coherent potential approximation is evaluated for comparision. In 
Sec. IV we discuss the concentration dependence of the critical temperature comparing 
predictions of the Belitz theory, the McMillan formula and the BCS theory. To this end 
we first establish density independent material parameters with the help of which one 
can construct the density dependent coupling strength A, the Coulomb pseudopotential 
/i and the Belitz correction Y' . Section V brings conclusions. 



2. Belitz Theory 



Let us outline the Belitz theory, first. In the absence of currents, the Eliashberg 
selfenergy has three terms E(e, u;) = (1 — Z(e, uj))ujtq + Y(e,u) f 3 + <f)(e,uj) f\ , where 
e is the kinetic energy and u is the Matsubara frequency. For pure superconductors </> 
renormalised by Z gives a gap function A while Y is a scalar that only shifts the chemical 
potential having no effect on the critical temperature. In disordered superconductors, 
parameter Y becomes energy and frequency dependent. Belitz has shown that its energy 
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is the screened pseudopotential, in which fi characterises a strength of the Coulomb 
interaction, ujq is its effective range and co>d is the Debye frequency. Setting Y' — one 
recovers the McMillan formula. 

In general, Belitz formula ([T]) includes four disorder dependent material parameters 
ub, A, fi and Y'. Based on ab-initio computations [HI [10] which show that oj^ changes 
negligibly with doping we ignore the disorder dependence of uj-q and use a pure diamond 
value hw-o/kB = 1860K. The strong disorder dependence of remaining parameters stems 
from a rapid variation of the density of states at the Fermi energy level N Q in vicinity 
of the metal-insulator transition. 
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We will use a simplified version of the Belitz theory neglecting disorder corrections 
to interaction vertices in perturbation series. This brings a significant reduction of 
numerical demands, because A, /i and Y' reduce to functions of Nq only. Moreover, with 
the vertex correction neglected, 

Y , = 2X~/2. (3) 

so that we are left with A and [i similarly as in the McMillan theory. 

We take A from ab-initio computations. Since the published results cover only few 
concentrations, we evaluate A in the spirit of Morel and Anderson [T2] formula 

T+~QN~o ' ( ) 

Here U represents a phonon-electron coupling strength and 1 + QNq describes a 
screening. 

The pseudopotential \i reads [U [13] 

where V is the strength of the Coulomb interaction. Since we have found no 
measurement of V in the literature, we have to treat it as a parameter which is set 
from experimental data. Because in the boron-doped diamond the Fermi energy level 
is localised very close to the top of the valence band we associate E-p « ujq with the 
effective range of the coulombic interaction and estimate both as the distance of the 
Fermi level from the top edge of the impurity band. The equation (j5J) then has a simple 
form 

fi = VN . (6) 
Using formulas (j3J), (BJ and (J6j) in (CD) and (j2J), we obtain T c as a function of x. 



3. Density of States 

To proceed we have to establish N Q as a function of the boron concentration x. This 
single particle property is independent of the electron-electron interaction, would it 
be the Coulomb repulsion or the phonon-mediated pairing potential. The density of 
states is thus given by a Hamiltonian of the valence band in diamond H and a random 
potential of boron impurities 

H = H + ^2ru5a%. (7) 

i 

In the random potential rji = 1 at impurity sites and zero elsewhere, and 5 is the 
potential amplitude. We note that this single-site s-type potential does not cover the 
triple degeneracy of impurity state of boron. We adopt this Hamiltonian already studied 
by Shirakawa et al [9] to make a link between their and present results. 

We have evaluated Nq within the DCA for clusters of 2 x 2 x 2 atoms on a cubic 
lattice. Compared to the CPA density of states employed by Shirakawa et al , our cluster 
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Figure 2. Density of states at the Fermi energy N as a function of boron 
concentration x computed using the CPA method and the DCA on the cluster 
of 2 x 2 x 2 atoms. The figure includes dashed curves N^ PA w 0.59 a; - 568 x 
10 23 cm- 3 cV _1 approximating CPA and Ng CA « 0.59 a; - 523 x 10 23 c m - 3 eV _1 for 
DCA. In most regions approximations are indistinguishable from computed results 
within the lincwidths. 
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density of states includes nontrivial corrections given by bonding and anti-bonding states 
at neighbour impurity sites. Difference between the CPA and DCA density of states in 
the impurity band is ilustrated in figure Q] for 5% of boron doping. 

Within the CPA (dashed line) the impurity band yields rather featureless density of 
states having a slightly skewed semieliptic shape. The valence band starts approximately 
at an energy of —1.5 eV and ends at —23 eV. We focus on the impurity band because 
the Fermi energy, E — 0, lies there. Since boron is an acceptor, its bound states have 
higher energy than extended states in a valence band. 

In the DCA (solid line) one can distinguish additional contributions of two and three 
boron clusters. Clusters of four to eight borons are also included, but their contributions 
are invisible for the given concentration. The main part of the impurity band density of 
states is formed by a bound state on a single boron. A shoulder on the right side of the 
impurity band results from a symmetric bound state of two neighbouring borons. The 
nonsymmetric state is not bounded but forms a resonant state visible in the valence 
band. The tree boron bound state is splited off at an energy 2.8 eV. Higher states of 
the three boron clusters make a negligible contribution. 

With respect to superconductivity the Fermi energy region is the most important. 
As one can see the density of states at E = is dominated by the single boron states. 
This is a reason why the CPA and DCA give comparable iVo- 

The resulting density of states at the Fermi level No is presented in figure [2] 
For comparison we also show the CPA result. One can see that the DCA density 
of states is higher than the CPA value at all impurity concentrations. Both resulting 
densities with a good accuracy obey power-law, N DCA « 0.59 x 0,568 x 10 23 cm _3 eV _1 
and Nq PA « 0.59 x ' 523 x 10 23 cm _3 eV _1 . A very good approximation by power- 
law also holds for the distance of the Fermi level from the top of the impurity band, 
LOc ~ 45.5 iVo 1 ' 18 eV for the DCA. In the rest of this section we provide more details of 
our computation. The reader not interested in these technical details can skip to the 
next section. 

We only briefly introduce the DCA here, the reader can find all details in 
reference [H]. A basic idea of the dynamical cluster approximation is to divide the 
Brillouin zone into subzones - in our case 2x2x2 subzones. Within the subzone the 
selfenergy is momentum-independent, i.e., the selfenergy is represented by eight complex 
functions of frequency. By symmetry arguments this number can be reduced to four 
functions. 

All subzones contribute to the density of states 

p°(e) = J2p k ( e ^ ( 8 ) 

K 

where K is a subzone index. In general, the subzone contribution p^ is obtained 
integrating over the subzone with the electron dispersion of the valence band. For 
simplicity we approximate these functions by a semielliptical distribution with a proper 
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Figure 3. Electron-phonon coupling A as a function of boron concentration x. Crosses 
are results of ab-initio computations using super cell method [HI [191 013 ■ Plus- 
signs represent virtual crystal computations [T51 [TTJ [T7] . Solid line is Morel- Anderson 
formula (gj) with U = 42.31 x 1CT 23 cm 3 eV and Q = 54.65 x 1CT 23 cm 3 eV. Densities 
of states have been associated to individual computations via DCA results shown in 
figure O 



width, shift and normalisation 

p>) = ^ ' /l-(^K) 2 , (9) 

where N c is number of subzones (8 in our case), and 

u K = l/2(£g- ax - £ff ) , (10) 

v K = 1/2(E^ + E™ 1 ) (11) 

are width and center of subzone energy bands with E^ x and E^ 1 being the maximum 
and minimum energy in the subzone K, respectively. 

The density of states (jHJ) by definition maintains the width of the valence band. 
Since one of the subzones is centered around the T point of the Brillouine zone, this 
approximation also yields a correct curvature at the edge, i.e., it correctly reproduces an 
effective mass of holes near the band edge. This feature is vital for a realistic description 
of the relatively shallow impurity state. We employ the model valence band (jHJ) in the 
CPA as well. 

The parameters of the hamiltonian (JTj) are fitted according to the facts that the 
valence band is 22 eV wide and the single impurity bound state appears at energy 0.37 eV 
above the valence one. Using the local Green function corresponding the density of states 

we come to 5 = 8.91 eV. 

We determine the Fermi energy level from a condition for a local averaged density 
of electrons at zero temperature n a = J°° p(E) dE = x/2. We are aware that a 
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Figure 4. Coulomb pseudopotential fi as a function of density of states N of boron 
doped diamond. Full (empty) symbols are values of fi deduced by Belitz (McMillan) 
formula from experimental data of Ref. [2J O , Ref. [2DJ 0, Ref. [T] V, Ref. [2T] A 
and Ref. [TS] □■ We employ A given by solid line in figure [3] The full symbols 
fall close to the solid line of the theoretical pseudopotential (JDJ with the Coulomb 
strength V — 2.54 x 10 - 23 cm 3 eV and the empty symbols are well reproduced with 
V = 8.18 x l(r 23 cm 3 eV. 

measured number of charge carriers can differ significantly from a number deduced 
from the concentration of boron atoms in the sample. As has recently been pointed out 
by Mukuda et al [15] quite a large fraction of boron atoms can form the neutral B-H 
complexes, reducing the concentration of charge carriers. On the other hand Klein et al 
[2] found that the effective number of carriers deduced from Hall-effect measurements 
was much larger than the number of boron atoms in samples. Since the charge carrier 
concentrations for all samples are not accessible we assume for simplicity the films to 
be doped ideally. 

4. Critical Temperature 

In this section we discuss the concentration dependence of the critical temperature T c . 
To be able to implement the Belitz theory we need the coupling strength A and the 
pseudopotential \i as functions of N . 

The coupling strength A we deduce from ab-initio computations. In literature one 
finds the studies within the virtual crystal approximation [161 EH D2] and the supercell 
method [HI [TU [TO]- In figure [3] we show fit of ab-initio results by formula @. We 
found that parameters U = 42.31 x 10~ 23 cm 3 eV and Q = 54.65 x 10 _23 cm 3 eV yield a 
reasonable fit of rather scattered computed values. 

It remains to establish the pseudopotential. According to formula ([6]) we have to 
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Figure 5. Critical temperature T c as a function of boron concentration x. Symbols 
are data as in figured] Thick (thin) solid line is a result given by the Belitz theory 
using the DCA (CPA) density of states. Results of the McMillan theory are in dashed 
lines. 

find a single parameter V, which determines \x for all concentrations. Dots in figure H] 
show values of \i deduced from the Belitz theory and experimental values of T c . One can 
see that all dots stay close to a line given by formula ([6]) with V = 2.54 x 10~ 23 cm 3 eV. 
The parameter V holds for all concentrations and it is the only material parameter fitted 
to experimental values of T c in the present theory. 

Now all relations and parameters are ready for predictions of the critical 
temperature. The concentration dependence of T c given by the Belitz formula (fTj) 
and power-law approximation of the DCA value of A*o is shown in figure El As one 
can see, the theory describes a steep increase of T c with doping at the region of small 
concentrations. At higher concentrations the critical tempeture satures at about 3 K. 
This saturation reminds trends described by the theory of Osofsky et al [22J based on a 
heuristic rescaling of the BCS parameters. 

For comparison we also show the value predicted by the Belitz theory for the CPA 
density of states. As one expects from the lower density of states, the CPA leads to lower 
values of T c , namely in the region of small concentrations. Our results thus confirm a 
trend predicted by Shirakawa et al that corrections beyond the CPA will lead to higher 
T 

J - c- 

We would like to emphasis that the disorder corrections by Belitz are necessary for 
a good agreement between experimental data and theory. It can be seen from attempt 
to the fit T c with the McMillan formula as shown in figure [5j In the comparison we have 
not used the pseudopotential /i fitted within the Belitz theory, but made a seprate fit 
directly from the McMillan formula, see figure HI 
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Figure 6. Failure of the BCS approach to the boron doped diamond. While 
experimental data (symbols) show decreasing slope of T c as a function of boron 
concentration x, the BCS theory (dashed lines) predicts steeply increasing slope. The 
parameter Vbcs = 3.97 x 10 _23 cm 3 eV was fitted within the DCA No computations. 
Thick (thin) line is given by the DCA (CPA) N Q results. Comparing the DCA (thick) 
and CPA (thin) one can see that the BCS theory overestimates effects of corrections 
beyond the CPA on the critical temperature. 



Finally we show that the boron-doped diamond has to be treated as a strong 
coupling superconductor. In figure [6] we show 100 experimental T c data and its 
description with the BCS formula T c = 1.14 exp (—1/Nq Vbcs) , where Vbcs = 3.97 x 
10 -23 cm 3 eV is the BCS interaction. Apparently, the BCS theory yields an incorrect 
concentration dependence. It is noteworthy how much the corrections beyond the CPA 
increase T c within the BCS theory. In the strong coupling theory these corrections are 
smaller by an order of magnitude. 

5. Conclusion 

To conclude we have discussed corrections beyond the CPA to the critical temperature 
in diamond doped with boron. From numerical studies of the dynamical cluster 
approximation with cluster of 2 x 2 x 2 atoms it follows that these corrections increase 
the density of states at Fermi level. As expected this causes an increase of the critical 
temperature. Comparing different approximations we have shown that the strong 
coupling theory is necessary to predict realistic critical temperature. It also turns out 
that disorder corrections of Belitz improve agreement with experimental data. 

The present theory offers a simple analytic formula for the critical temperature as a 
function of boron density. Indeed, in a simplified Belitz theory all material parameters 
depend exclusively on the density of states and can be expressed as its linear or rational 
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functions. Numerical studies of C\- X B X show that the density of states at the Fermi 
level follows a power-law N ~ 0.591 x a523 x 10 23 cm~ 3 eV _1 , what allows one to predict 
the critical temperature without any numerical effort. 
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